J/ip polarization at hadron colliders in nonrelativistic QCD 
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With nonrelativistic QCD factorization, we present a full next-to-leading order computation of the 
polarization observable for J/ip production at hadron colliders including all important Fock states, 
i.e. ^jj 1,8 ' , ^Sq 8 ' , and ^j 8 ' . We find the ^j 8 ' channel contributes a positive longitudinal component and 
a negative transverse component, so the J/ip polarization puzzle may be understood as the transverse 
components canceling between %*[ 8 ' and ^j 8 ' channels, which results in mainly the unpolarized (even 
slightly longitudinally polarized) J /ip. This may give a possible solution to the long-standing J/ip 
polarization puzzle. Predictions for J/ip polarization at the LHC are also presented. 
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Nonrelativistic QCD (NRQCD) @] is an effective field 
theory approach for heavy quarkonium. At present, one 
of the main obstacles to NRQCD is the polarization puz- 
zle of J/ip hadroproductionQ. At leading order (LO) in 
a s , J/ip production is dominated by gluon fragmentation 

to a color-octet (CO) S[ cc pair at high transverse mo- 
mentum pt , which leads to transversely polarized J/ip || . 
But the CDF Collaboration found the prompt J/ip in its 
helicity frame to be unpolarized and even slightly lon- 
gitudinally polarized H). Despite of numerous attempts 
made in the past, the puzzle still remains. 

For unpolarized J/ip production, important progress 
has been made in recent years. It was found that the 
next-to-leading order (NLO) QCD corrections to diffcr- 
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cntial cross sections of S[ channel can be as large as 2 
orders of magnitude at high pr H , while those of the ^Sq 8 ' 
and sj 8 ' channels are small Q. Furthermore, NLO cor- 
rections of the 7j [jjj and 3 pf^ || [| channels are found 
to also be very large. These large corrections are well 
understood because at NLO the differential cross sec- 
tion dcr/dpT receives contributions from new topologies 
that scale with px in a different manner from the LO 
calculation. By including NLO corrections, one may ex- 
plain the existing unpolarized cross sections of px up to 
70 GeVEJ. 

Accordingly, it is necessary to examine the J / ip polar- 
ization at NLO. Among various channels, the correction 
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to J/i p polarization via the S[ channel was worked out 
in |lCf| , which alters the polarization from being trans- 
verse at LO to longitudinal at NLO. This phenomenon 
was explained recently in collinear factorization pd| as 
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next-to-lcading power dominance. As for the S[ chan- 
nel, the NLO correction can only slightly change the 
polarization^, while the Sq 8 ' channel gives an unpolar- 



ized result to all orders in a s . As a result, up to date 
theoretical predictions may indicate a serious puzzle for 
J/ip polarization[||. However, the NLO correction of the 
^-j 8 ' channel to J / ip polarization has not been calculated. 
In this Letter, we perform this calculation and show that 
the NLO contribution of channel is indeed crucial in 
clarifying the long-standing J/ip polarization puzzle in 
NRQCD. 

We first introduce some formalisms in our calculation. 
The J/ip can decay into an easily identified lepton pair. 
The information about the J/ip polarization is encoded 
in the angular distributions of the leptons. The two- 
body leptonic decay angular distribution of J/ip in its 
rest frame is usually parameterized as]l2| 
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Here, day = 0, ±1, with respect to the z compo- 
nents of J/ip) represents the ij contribution in the spin 
density matrix formalism. In the literature, the polariza- 
tion observable \g is also denoted as a = ^t-Mjl -phe 
differential cross sections are 

dcr s zSz =][] j dxidx 2 / i/Hi (xi, pi F )/ j/H2 (x 2 , /i F )(O n )d^ J ; s n a ,(2) 

where (On) are the long-distance matrix elements 
(LDMEs) for n = 3 sj 1,81 , ^j 81 , and !sf ] . In general, the 
partonic cross sections da^ can be obtained from the 
spin density matrix elements fl2| 

P.„..(*j -> (cc)[n]X) cx ^\M(ij -> (cc)[L z ,s z ]X)\ 2 .(3) 



In practice, several polarization frame definitions have 
been used in the literature. In the s-channel helicity 



2 



frame, the polar axis is chosen as the flight direction of 
the J/ip in the laboratory frame. Another frequently 
used frame is the so-called Collins-Soper frame [|l3]. For 
simplicity, here we will only choose the helicity frame, the 
same as used by CDFQ]. The full theoretical predictions 
of azimuthal correlations and the theoretical descriptions 
by Collins-Soper, and feed down from \cj and ip' will be 
presented in a forthcoming publication. 

We describe our method briefly for the sake of com- 
pleteness. Some improvements are made in our calcula- 
tions, while most of our method has been encompassed 
in Ref . [B . The calculations of real corrections are based 
on the Dyson-Schwinger equations. After absorbing the 
core codes of the published HELAC|l4|], we promote it 
into a form that can generate the matrix element of 
heavy quarkonia (especially P-wave) production at col- 
liders by adding some P-wave off-shell currents. The vir- 
tual corrections are treated analytically, and the helicity 
matrix elements are obtained using the spinor helicity 
method @. 

For numerical results, we choose the same input pa- 
rameters as in Ref. [||. Specifically, the renormalization 
scale fi r , factorization scales [if and NRQCD scale /it a are 
chosen as fi r = /if = mx = \J 4m 2 + p\ and fj>\ = m c . 
Scale dependence is estimated by varying fi r ,fif, by a 
factor of ^ to 2 with respect to their central values. By 
fitting only the cross sections, it was found that only two 
linear combinations of CO LDMEs can be cxtractcd[||. 
Since polarization information is also available, we will 
try to extract the three independent CO LDMEs using 
the polarization observable Xg and the production rate 
dcr/dpT of the J ftp measur cd by CDF Run II g simulta- 
neously, where the data in the low transverse momentum 
region (px < 7GeV) are not included in our fit because of 
existing nonperturbative effects. By minimizing x 2 , the 
CO LDMEs are obtained and listed in the first row of 
Table.0. In Fig.ffl we compare Xg from the Tevatron data 
with our theoretical results. 



GeV 3 


10 _2 GeV 3 


(C(^l 81 )} 
l(T 2 GeV 3 


10 _2 GeV 3 


1.16 


8.9 ±0.98 


0.30 ±0.12 


0.56 ±0.21 


1.16 





1.4 


2.4 


1.16 


11 









TABLE I: Different sets of CO LDMEs for the J/ip. Values 
in the first row are obtained by fitting the differential cross 
section and polarization of prompt J/ip simultaneously at the 
Tevatron Q . Values in the second and third rows are two ex- 
treme choices for these CO LDMEs. The color-singlet LDME 
is calculated by the B-T potential model in |Tt| ). 



To understand the unpolarizcd results, Xg for each 
channel is drawn in Fig.|^, where for the NLO 3 pj 8 ' chan- 
nel we mean the value of (d<xn — d<7oo)/|d<7n + d<7oo| 
because d&n + dcroo decreases from being positive to 
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FIG. 1: (color online) NLO results for the polarization ob- 
servable Xe of J/ip production at the Tevatron. The CDF 
experimental data are taken from Ref. W. 



negative as px increases. In addition to the known 
polarization of S-wave|], [[(J, the 3 pf^ channel satisfies 
(dan — d(7oo)/|d(Tii ± doool < — 1 in our considered px 
region, which results from dan < and d<7oo > 0. There- 
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fore, the transverse component of Pj is negative, which 
effectively gives a longitudinal contribution to Xg, and 

3 roi 

the longitudinal component of Pj is positive. In some 
parameter space of the CO LDMEs, the positive trans- 
verse component of S[ will largely be canceled by the 

3 T8l 

negative transverse component of Pj , which yields a 
small transverse component and results in an unpolarized 
or even longitudinal Xg. This explains why the complete 
NLO calculation gives an unpolarized prediction in FigJ^. 

It is interesting to see that, by choosing the proper 
CO LDMEs, complete NLO predictions in NRQCD fac- 
torization can be made compatible with the data. This is 
distinct from all previous NRQCD predictions that give 
strong transverse polarizations for J/^[g. Furthermore, 
we want to emphasize the following four points. 

(1) Transverse components with large cancellation be- 
tween and Tj 8 ^ determine a specific parameter space 
for the CO LDMEs. Using the same treatment as in 
Ref. Q|, we decompose the transverse component of the 
short-distance coefficient of 3 p| 8 ' into a linear combina- 
tion of 3 Sf and as dan^Pf) = 2.47 dan^ 81 ) - 

0.52 danOsf 1 ). Since dtrn^ 81 ) < dan^f 1 ) when 
Pt > 7 GeV, the cancellation requirement is approx- 
imately equivalent to the absence of the linear combi- 
nation (0( 3 Sf ] )) - 0.52 (0( 3 P [ f ] ))/m 2 c , which is close to 

Mr = (0(M 81 )) - 0-56(C( 3 p{ r 8l ))/m2 defined in Ref.f). 
Recall that to have a good fit for the unpolarized yield 
one needs a very small Mi, so the conditions for the CO 
LDMEs parameter space introduced by fitting both yield 
and polarization are consistent with each other. Good 
agreement with the LHC data for the J/ip cross sections 
can be found in Fig|3| using the LDMEs in Table|j|. 
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(2) As the yield and polarization share a common pa- 
rameter space, and the yield can only constrain two linear 
combinations of CO LDMEs, the combined fit of both 
yield and polarization may also not constrain three in- 
dependent CO LDMEs stringently. In fact we find for 
a wide range of given (O(S l ')), one can fit both yield 
and polarization reasonably well. CO LDMEs under two 
extreme conditions are listed in Tablef When (C>Cs [ * ] )) 
is chosen to be its maximal value, the J/tp is unpolar- 
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ized; when (0( Sq )) vanishes, Xg increases from -0.25 at 
Pt=5 GcV to at pt=15 GeV at the Tevatron. Even in 
these two extreme cases, the theoretical predictions of the 
J/tp cross section and polarization are still close to the 
Tevatron data, and are also consistent with the observed 
cross sections obtained by ATLAS @ and CMS@ at 
the LHC as shown in Fig. | As a result, although it is 
hard to determine the CO LDMEs precisely, we find that 
the polarization puzzle can be much eased for a wide 
range of (0(^j 8] )} value. 

(3) The cancellation of the transverse component be- 
tween the 3 S^ and ±j channels is not problematic, 
since the contribution of an individual channel is un- 
physical and depends on rcnormalization scheme and 
scale ||. A "physical" requirement is that the summa- 
tion d(7ii('fe!j 8 ' + "Pj 8 ') be positive, which is satisfied in 
the fit. 
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FIG. 2: (color online) The pr dependence of Xg for 

^f 1 ,^ 81 ,^ 81 , and ^j 81 channels with VS = 1.96 TeV and 

\yj/i>\ < 0.6. For the NLO Vf^ channel, the dot-dashed curve 
means the value of (d<rn — d<7oo)/|d<7ii + daoo|- 



(4) It is important to note that the LDMEs presented 
here are significantly different from those extracted from 
the global fit in Ref. . As hadroproduction data play 
the most important role in Ref.Jl|], this difference can- 
not be mainly attributed to data other than hadropro- 
duction not being considered in our fit. In fact, one 
can track to the situation where only hadroproduction 
data are used in the global fit. As explained in Ref.(8|, 
our choice of the pt cutoff for hadroproduction data is 



Pt > 7GcV while the cutoff in 1 19 is px > 3GcV, and our 
LDMEs can well describe the pt spectrum in the region 
7GcV < p T < 70GcV (sec Fig.|), while the fit in Ref.Q 
puts stress on the smaller px region and gives too smooth 
a pt distribution at large pt ■ This is the main reason why 
our LDMEs differ from those in Ref. |l9| . In our view, for 
the small px region the fixed order perturbation calcula- 
tion may need to be modified by considering soft gluon 
emission and other nonperturbativc effects. We see that 
the two treatments in Ref. || and Ref.jl9| have different 
features and should be tested by more experiments in the 
future. 



and the central region (\yj/^\ < 2.4) arc considered 1 . 
The large error bar (yellow bands) in these predictions is 
caused by a lack of knowledge of (0( S )); thus, we scan 
all its possible values in the predictions. It is found in 

1 f8l 

these predictions that Xg becomes sensitive to (0(Sq )) 
whenpr > 20 GeV, so it may be possible to extract three 
independent CO LDMEs when polarization data at high 
Pt are available. 

In summary, we present a full NLO calculation includ- 
ing sj 1 ' , 3 S^ , "feg 8 ' and 3 pj 8 ' for the polarization observ- 
able Xg of the J I ip in the helicity frame at the Tevatron 
and LHC. Results of S-wave channels are consistent with 
those in the literature |l(| , while those of the ^Pf^ channel 
are new and play a crucial role in understanding the po- 
larization puzzle. Our calculation shows that the trans- 
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verse component of the Pj channel is negative, while 
its longitudinal component is positive. Thus the 

3p[8] 

channel gives a maximal longitudinal contribution. By 
choosing suitable CO LDMEs, which bring about good 
agreement with the observed J/tp cross sections at large 



There are still other uncertainties, such as the charm 
quark mass, but they do not change the qualitative prop- 
erties of our result. Predictions of the polarization ob- 
servable Xg at the LHC with \/~S = 7 TeV are plotted 
in Fig.^J, where only the forward region (2 < I2/J/0I < 3) 



1 Note that the ALICE Collaboration has measured J/ip polariza- 
tion recently with rapidity 2.5 < \yj/^\ < 4[E(J. But the mea- 
sured transverse momenta (2GeV < pt < 8GeV) are smaller 
than considered in this Letter. 
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FIG. 3: (color online) Cross sections of J/ip production at the 
LHC with LDMEs shown in Table |. The green bands (NLO 
Total I) correspond to the LDMEs in the first row of Table.|. 
The yellow bands (NLO Total II) correspond to the LDMEs 
in the second row 

[(C( 1 4 8] )) = (upper bounds)] and third 
[(ffC&l 81 )) - (^(^Pj 81 )) ~ ( lower bounds)] of Table 



Pt at the LHC, the transverse components can be largely 
canceled between the 3 S^ and ^j 8 ' channels, leaving 
the remaining terms to be dominated by the unpolar- 
ized J/i/). This may give a possible solution to the long- 
standing J/tp polarization puzzle within NRQCD factor- 
ization. Although it is hard to individually extract the 
three independent CO LDMEs in an accurate way, our 
interpretation of J/)j: polarization makes sense by using 
only their combinations. We also present polarization 
predictions for the LHC. 
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Note added. When this Letter was being prepared, a 

preprint [ BjJ on the same issue had just appeared. The 

— 3 rsi 

essential difference is that they have a negative (0( Pq )) 

based on a global fit , and give a significant transverse 
polarization prediction, but our fit leads to a positive 
(0(Pq )), which is consistent with observed cross sec- 
tions in a wide pr region (7-70 GeV) at the LHC and 
results in the mainly unpolarized J/ip. 
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FIG. 4: (color online) NLO predictions of the J /if) polariza- 
tion observable Xe at the LHC. The uncertainty is shown by 
large yellow bands when varying the CO LDME (C(^ 81 )). 
The bounds of \q — in yellow bands correspond to CO 
LDMEs in the third row of Table.|, while the other bounds 
correspond to the second row of Table. |. The small green 
bands are the predictions using the CO LDMEs in the first 
row of Table|. 
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